Conceptualization of seeded region growing by 
pixels aggregation. Part 3: a wide range of 

algorithms. 
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Abstract — In the two previous papers of this serie, we have 
created a library, called Population, dedicated to seeded region 
growing by pixels aggregation and we have proposed different 
growing processes to get a partition with or without a boundary 
region to divide the other regions or to get a partition invariant 
about the seeded region initialisation order. Using this work, 
we implement some algorithms belonging to the field of SRGPA 
using this library and these growing processes. 

Index Terms — Distance function, dynamic filter, geodesic re- 
construction, homotopic transformation, regional minima, seeded 
region growing by pixel aggregation, vorono'i tessellation, water- 
shed transformation. 



I. Introduction 

Many fields in computer science, stereovision[12], matli- 
ematical morpiiology[14], use algoritiim wiiicii principle is 
Seeded Region Growing by Pixels Aggregation (SRGPA). 
This method consists in initializing each region with a seed, 
then processing pixels aggregation on regions, iterating this 
aggregation until getting a nilpotence [1][10]. The general 
purpose of this field is to define a metric divided into two 
distinct categories [3]: the region feature hke the tint [1] and 
region boundary discontinuity [6]. 

In this article, the aim is not to do an overview of the 
algorithms using SRGPA but to prove that the framework 
introduced in the two previous articles [15] [16] is generic. 
Some algorithms using SRGPA are implemented thanks to the 
library Population; 

• voronoi tessellation, regional minima, domain to clusters, 

• distance function, watershed transformation and geodesic 
reconstruction. 

The first enhancement is the easiness to implement these 
algorithms using the objects of the library Population. The 
second enhancement is the algorithms efficiency. All these 
algorithms have been applied on 3D image with a size equal 
to 700*700*700=0.35 Giga pixels. The running time is always 
less than 3 hours with an Intel(R) Xeon(R) CPU 3.00GH. This 
is due to 

1) the library optimisation using the template metaprogram- 
ming'[2]; all algorithms using this library will benefit 
from this optimization, 

'Template metaprogramming is a metaprogramming technique in wliicli 
templates are used by a compiler to generate temporary source code, which is 
merged by the compiler with the rest of the source code and then compiled. 
The output of these templates include compile-time constants, data structures, 
and complete functions. The use of templates can be thought of as compile- 
time execution. 



2) the procedure of actualization of the zones of influence 
described in the previous article[15]. 
In this article, the notations are: 

• let E' be a discrete space^, 

• let O be a domain of E and / its characteristic function 
such as n ^ {\fx G E : I{x) ^ 0}, 

> let / be a grey-level image, an application of E to Z, 

> let F be a neighborhood function (an elementary struc- 
turing element). 

In the appendice I, the definition of the distance is given. 
The article understanding depends on the comprehension of 
the previous articles of this serie. A summary is done in the 
appendice II. 

The outline of the rest of the paper is as follows: in Sec. II, we 
present the algorithms using only one queue in the system of 
queue (SQ), in Sec. Ill we present the algorithms using more 
than one queue, in Sec. IV, we make concluding remarks. 

II. One queue 

In this section, we will present some algorithms using a 
single queue during the growing process. 



A. Simulated Voronoi tessellation 

Consider $ a Poisson point process in a metric space M. 
The cells 

C{x) ^ {y e M; d{y ~ x) < d{y - x), x e $}, a; e $, 

constitute the so-called Poisson- Voronoi tessellation of M. 
Presented by Gilbert in 1962 [8], this statistical model is 
appropriate for random crystal growth. In the discrete space 
E, the implementation for a distance associated to norm 1 or 
cx)^[13] is done using the library Population. 
Starting form the affectation of each region with a seed (a 
point of Poisson point process), an isotopic growing process 
at constant velocity is operated. The ordering attribute function 
is b(x, i) = 0. The growing process is (see algorithm 1 and 
figure 1): 

> initialization of the regions/ZI by the seeds 

• select the queue number 

* while the selected queue is not empty 

^The space E, is a n-dimensional discrete space Z", consisting of lattice 
points whose coordinates are all integers in a three-dimensional Euclidean 
space R". The elements of a n-dimensional image array are called points. 

'For the Euclidian distance, see [17]. 
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- extract {y,i) from the selected queue 

- "Growth on x of the region i" 

• return regions 

The quote "growth on x of the region i" means that there 
are different kinds of growing process introduced in the 
previous article [16]. Here, the growing process is done 
without a boundary region to divide the other regions. In 
the algorithm 1, the growing process leading to a final 
partition invariant about the seeded region initialisation order 
is used. To prove that this growing process gives a correct 
Poisson-Voronoi tessellation of E, this property is used: 

Va;, y G E : d{x, y) = min {(d{x, z) + 1) 

The generation of a Poisson point process is done using the 
Boost software. This implementation is not restricted to the 
Poisson-Voronoi tessellation since: 

• each seed can be a domain of E (second serie in the 
figure 1), 

• the growing process can be restricted to a domain O = 
{Vi G E : I{x) 7^ 0} if the ordering attribute function 
is: 5{x, i) = if I{x) ^ 0, OUT else (third serie in the 
figure 1). 



Algorithm 1 Geodesic dilatation with an invariant boundary 
Require: S , V //The binary image, the seeds, die neighbor- 
hood 

II initialization 

System.Queue s_q( 6{x, i) 0, FIFO, 1); //A single FIFO 
queue 

Population p (s_q); //create the object Population 

Tribe passive(y = 0); 

int ref_boundary = p.growth_tribe(passive); 

//create a boundary region/ZI, (X^, Z^) such as Z\ = 

Restricted A^=N; 

Tribe active(V, N); 

. Vsi e S" do 

int ref_tr = p.growth_tribe(actif); //create a region/ZI, 

{XI ZD such as Z\ = {X\ ®V)\{\} X,) 

J en 

p.growth(si, ref_tr ); 
end for 

//the growing process 

s_q.select_queue(0); //select the single FIFO queue. 
while s_q.empty()==false do 
{x, i) = s_g.pop(); 
if pop.Z()[x].size()>2 then 

p.growth(x, re f .boundary); //growth of the boundary 
region 
else 

20: p.growth(x, i ); //simple growth 
end if 
end while 
return p.X(); 
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Fig. 1. Each serie serie is the Voronoi-growing process. For the first 
serie, the seeds are drawn from the Poisson point, A video is avail- 
able at http://pmc.polytechnique.fr/^vta/geodesicJmariant_cube.mpeg. For 
the second serie, each seed is a set of points, A video is available at 
http://pmc.polyteclmiqiie.fr/~vta/geodesic. imariant.cube.grain. mpeg. For the 
last serie, the growing process is restricted by the first figure, A video is 
available at http://pmc.polytechnique.fr/~vta/geodesic.mpeg. 



B. Domain to clusters 

Let O be a domain of E and let be the set of continuous 
application from [0, 1] to VI. 
(ci)o<(<n is the clusters decomposition of if: 

Uo<i<nCi = Vl 

e (0,...,n)Vx,2/ g (c,,q) 37 gC" 7(0) = x A 7(1) = 2/ 
Vz^jVa;,y e (q,Cj) ^ ^ 7(0) = x A 7(1) = y 

The second line means that all points belonging to the same 
connected component are linked and the third line means 
that two points belonging to different connected components 
are not linked. This extraction gives information about the 
critical percolation concentration, percolation probabilities, 
and cluster size distributions[ll]. Using the library Popula- 
tion, a algorithm is defined to extract the set of connected 
components. The principle is: when a connected component 
is touched, this connected component is removed from Q, using 
a growing process (see algorithm 2 and figure 2): 

• scan the image (Vx G E) 

- if I{x) ^ 

* create a region/ZI initialised by the seed {a;} 

* select the queue number 

* while the selected queue is not empty 

• extract (y, i) from the selected queue 

• growth of the region i on y 

■ Hy) = 

• return regions 

Using this extraction, it is possible (see figure 3): 

• to remove all the connected components touching the 
boundary, 

. to fill the hole"*, 

• to keep only the cluster which area is maximum, 

"^To file the hole, the porcedure is 

1) inversion of the initial image, 

2) extraction of connected components, 

3) removing the connected components no touching the image boundary, 

4) binarization and inversion of this last image. 



Fig. 2. Left image: scan the image until I(x) 7^ 0, middle image: growing 
region starting form x (I{x) 7^ 0) such as at each growth the characteristic 
function of Q is modified I{x) 7^ ^ I{x) = 0, right image: at the end of 
the growing region, the connected component has been extracted and removed 
from n and the scanning continues until I{x) 7^ 0. 



Algorithm 2 Domain to clusters 



Require: /, V //The binary image, the neighborhood 
II initialization 

System.Queue s_q( 6{x,i) = if /(x) ^ 0, OUT else, 

FIFO, 1); //A single FIFO queue such as if I{x) = then 

(a;, i) is not pushed in the SQ. 

Population p (s_q); //create the object Population 

Restricted N ==N; 

Tribe active(V, N); 

llScan the image 

.yxeE do 

//Test if a connex component is touched 

if I{x) ^ then 

int ref_tr = p.growth_tribe(actif); //create a region/ZI, 
{XI ZD such as Zl = {XI ®V)\{[jX,) 

p.growth(x, ref_tr ); 
lithe growing process 

s_q.select_queue(0); //select the single FIFO queue. 
while s_q.empty()==false do 

(y, i) = s_g.pop(); 

p.growth(?/, i ); 

l{v) = 0; 
end while 
end if 
20: end for 

return p.X(); 



C. Regional minima 

Let C5,, be the set of continuous application from [0 



-x.y 

to E such as the two extremities are equal to 

(V7 e C^^y : 7(0) = X A 7(1) = y). 

S = (si)o< j<Ti is the decomposition of {E, I) in level 
connected sets if: 



1] 

and y 



U s., 

0<i<n 



E 



Vi G (0, . . . , n)Vx, y & Si 
V* jVj e Cf e [0, 1] 



Fig. 3. Left upper image: the initial binary image, left middle upper image: 
extraction of the connex components, right middle upper image: the connected 
components touching the boundary of the image are removed, right upper 
image: the max cluster of the previous image, left bottom image: the initial 
image, right bottom image: maximum cluster of percolation after the selection 
of the component whose the area is maximum in the connected components. 



minimum if: 

W{x,y) e {s,{s(BV)\s) 



I{x) < I{y) 



I{x) = I{y) 
(J(x)^/(7(t)))V(/(2/)^J(7(i))) 



If / is seen as a topographic surface, the second line means 
that the level is the same in each point belonging to Si and the 
third line means that all paths between two points belonging 
to different elements of S do not have a constant level. 
In this decomposition, an element s of 5 is a regional 



The level of the points belonging to the outer boundary of 
s is greater than the level of the points belonging to s (see 
figure 4). Using the library Population, a growing procedure is 
defined to extract the regional minima. This growing procedure 
consists to scan the image (Vx G E). At each time, there is 
not yet a region on x {pop. X{)\x\. empty {) == true) to start 
the growing region initialized by the seed equal to {a;}. Let 
level = I{x) be the level of the growing region. The ordering 
attribute function is defined as: 

5{x, i) = if I{x) < level, OUT else 

For this algorithm, the ZI is defined as: Z* = (X- © V^) \ X* 
because the ZI is localized on the outter boundary region even 
if there are still some region to check the condition: V(a:, y) £ 
(s, (s © y) \ s) : I{x) < I{y). The growing process is (see 
algorithm 3 and figure 4): 

• to scan the image (ix G E) 

- if pop. X{)[x].empty{) == true 

* create a region/ZI initialised by the seed {x} 

* level = I(x) 

* select the queue number 

* while the selected queue is not empty 

• extract {y,i) from the selected queue 

• if I{y) == level 

■ then growth of the region i on y 

■ else 

• then this region/ZI is not a regional minimum 
return regions that are regional minima 



III. 



N QUEUES 

In this section, we will present some algorithms such as a 
muti-queue is used during the growing process. 
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Decomposition into level connected sets 
Legional minima 




Fig. 4. Upper image: principle of minima; bottom images: on tiie left, a 
grey-level image, on the right the regional minima of this image 



Algorithm 3 Regional minima 

Require: /, V //The grey-value image, the neighborhood 
II initialization 

int level; 

System.Queue s_q( 6{x, i) = if I{x) 
level, OUT else, FIFO, 1); //A single FIFO queue 
Population p (s_q); //c reate the object Population 



< 



Restricted N = {i}; 
Tribe active(V, N); 
Set set; //Container: self-balancing binary search tree. 

llScan the image 
.Vxe E do 

//Test if there is still a region on x 
if pop. X{)[x\.empty{) == true then 
level = I{x) 

int ref_tr = p.gro wth_tribe(actif); //create a region/ZI, 



{Xf, Zj) such as 
p.growth(x, ref_tr 
bool regional_minima=true; 
lithe growing process 



zl = {X\ © V) \ x\ 



X 



s_q.select_queue(0); //Select the single FIFO queue. 
while s_q.empty()==false do 

(y, i) = s_g.pop(); 

if I{y) < level then 
regional_minima=false; 

else 

p.growthCy, i ); 
end if 
end while 

if regional_minima==true then 

set.insert(ref_tr); 
end if 
end if 
end for 

return (p.X(),set); 



A. Distance function: flip-flop queue 

For the voronoi tessellation, we impose only a growing 
process at a constant velocity with forgetting the distance 
between the seeds and the current position. Here, it is a 



growing process step by step where a step corresponds to a 
distance value. The output is a distance function. 

a) n queue implementation: Let d be the distance be- 
tween a point and the seeds. The ordering attribute function 
is: 5{x, i) = d + 1 ^. The growing process is: 

. int d=0 

« initialization of the regions/ZI by the seeds 

• while the system of queues is not empty 

- d = d H- 1 

- select the queue number d 

- while the selected queue is not empty 

* extract (y, i) from the selected queue 

* growth on x of the region i 

* dist[x]=(i 

• return dist 

The number of queues is equal to the maximum of the distance 
function. The problem of this implementation is that this 
number is unknown before the growing process. To overcome 
this problem, a solution is to use a flip-flop queue. 

b) flip-flop queue implementation: In the last implemen- 
tation, during the growing process, there are only two queues 
in the SQ not empty at the step d: the queue number d 
where the couples are extracted and the number d + 1 where 
the couples are stored. Using this property, the couples are 
now extracted from the queue number flip and stored in the 
queue number flop. The ordering attribute function 5{x,i) is 
equal to flop. The growing process becomes (see figure 5 and 
algorithm 4): 

. int d=0 

• initialization of the regions/ZI by the seeds 

• while the system of queues is not empty 

- d = d H- 1 

- switch(flip,flop) 

- select the queue number flop 

- while the selected queue is not empty 

* extract (y, i) from the selected queue 

* growth on x of the region i 

* dist[a;]=(i 
« return dist 

This algorithm is not limited to the distance function of E. 
The growing process can be restricted to a domain Q, = {Vx e 
E : I{x) 7^ 0} if the ordering attribute function is: 5{x, i) = 
flop if I{x) ^ 0, OUT else (fourth serie in the figure 1). 

B. The watershed transformation 

An efficient segmentation procedure developed in mathe- 
matical morphology is the watershed segmentation [6], usually 
implemented by a flooding process from labels (seeds). 
Any greyscale image can be considered as a topographic 
surface and all boundaries as sharp variations of the grey 
level. When a gradient is applied to an image, boundaries 
are enhanced. When the topographic surface obtained from 

^The growing process is on the points which value is equal to d and each 
couple (x, i) is stored in the queue number d+ 1. 
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Fig. 5. The distance function 



Algorithm 4 Distance function 



Require: I, S , V //the binary image, the seed, the neighbor- 
hood 

II initialization 

Image Img_dist; 
int flip=0, flop =1; 

System_Queue s_q( S{x, i) = flip if I{x) ^ 

0, OUT else, FIFO, 2); //two FIFO queues. 

Population p (s_q); //create the object Population 

Restricted = N; 

. Vsi G S* in the order 0, 1 ... do 

int ref_tr = p.growth_tribe(actif); //create a region/ZI, 
[Xf, ZD such as Z\ = {X\ ®V)\{\]Xj) 

p.growth(si, ref_tr ); 
end for 
int dist=0; 

lithe growing process 

while s_q.alLempty()==false do 

switch(flip,flop); 

s _q. select_queue(flop) ; 

dist++; 

while s_q.empty()==false do 
{x, i)=s_q.pop(); 
p.growth(a;, i ); 
Img_dist(a;)= dist; 
end while 
end while 
return Img_dist; 



the gradient is flooded from its seeds, the waterfronts meet 
on watershed lines in 2D, and on watershed surfaces in 3D. 
A partition of the investigated volume is obtained, where the 
catchments basins are separated by the watershed surfaces (see 
figure 7). 

To implement this algorithm, a step by step growing process 
is defined where a step corresponds to a level of immersion. 
Let level be the level of immersion. The ordering attribute 
function is: S{x,i) = max(^eue^, f{x)) such as all the points 
immerged at the same level are stored in the same queue. The 
growing process is: 

• int level=f.min_range(); 

• initialization of the regions/ZI by the seeds. 




Fig. 6. First serie: the left image is a realisation of a random point process 
with A = 0.005 (for the visualisation convenience, the realisation has been 
dilated), the middle image is the regions at the end of the growing process, 
the right image is the distance function. Second serie: the same as the first one 
except the seeds are not some single points but some connected components. 
Third serie: the same as the first one except the seed is the complementary of 
the domain. Fourth serie: the first image is the domain and the second image 
is the distance function on this domain such as the seed has been localized 
on the blue face. This distance function is used to calculate the geometrical 
tortuosity. 



• for level = f.min_range() to f.max_range() 

- select the queue number level 

- while the selected queue is not empty 

* extract {y,i) from the selected queue 

* "growth on x of the region i" 

• return regions 

The quote "growth on x of the region i" means that there are 
different kinds of growing process introduced in the previous 
article [16]. In this example, the growing process is done 
without a boundary region to divide the other regions. In the 
algorithm 5, the growing process leading up to a final partition 
invariant about the seeded region initialisation order is used. 
This growing process is not limited to the watershed trans- 
formation on E. The growing process can be restricted to 
a domain = {Vx G E : I{x) ^ 0} if the ordering 
attribute function is: S{x,i) — ma,x{level, f{x)) if I{x) ^ 
0, OUT else (see figure 8). 
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Fig. 7. This transformation requires two images: the topographic surface 
(a grey-level image) and the label image (the seeds). The process is: 1) 
association of each label to a hole 2) immersion 3) the water enters in the 
topographic by the holes and the catchment basins take the colour of the hole, 
4) a part of the topographic is not merged although its level is under the level 
of the immersion, 5) fusion of two catchment basins with same colour, 6) 
creation of dam when two catchment basins have different colours. A video 
is available at http://pmc.polytechnique.fr/~vta/water.mpeg 




Fig. 8. The first serie: the first image is the initial image, the second image 
is the application of a Deriche's gradient[7] on the first image, the third image 
is the visuahzation of two seeds: one localized on the grains, the other on 
the grains complementary (in this case, the seeds are not connected). The 
second serie: the first image is the catchment bassins due to the watershed 
transformation on the previous gradient image using the two seeds, the second 
image is the foreground of the boundary region on the initial image. There is a 
good match with the visual segmentation. The third serie: the first image is the 
initial image, the second image is a slice of the opposite distance function of 
the initial image, the third image is the catchment bassins due to the watershed 
transformation restricted by the initial image on the opposite distance function 
of the initial image with appropiiate seeds. 



Algorithm 5 The watershed transformation with an invariant 
boundary 

Require: f,S,V //the topographic image, the seed, the 
neighbourhood 
II initialization 

int level=Q; 

System.Queue s_q( S{x,i) = uia.x{f{x),level), FIFO, 

/.max_range() - /.min_range()+l); //n FIFO queues. 

Population p (s_q); //create the object Population 

Restricted N = W, 

Tribe passive(y = 0); 

int ref_boundary = p.growth_tribe(passive); 

Tribe active(V, N); 

. Vsi e 5 in the order 0, 1 ... do 

int refJr = p.growth_tribe(actif); //create a region/ZI, 
ZD such as Z\ = (X\ ®V)\{\^X,) 

p.growth(si, refJr ); 
end for 

lithe growing process 

. For level = f.min_range() to f.max_range() do 

s_q.select_queue(^eweO; //Select the queue number level 
wliile s_q.empty()==false do 
{x, z)=s_q.pop(); 

if pop.Z()[x].size()>2 and i= min_elements( 
pop.Z()[x]) then 

p.growth(x, ref_boundary); //growth of the bound- 
ary region 
20: else 

p.growth(x, i ); //simple growth 
end if 
end while 
end for 

return pop.X(); 



C. Geodesic reconstruction 

The geodesic reconstruction is an efficient tool in Morphol- 
ogy Mathematic [14], [4]. Given a function / and a function 
g with / > g, the geodesic erosion is defined as: 

where E^{f) is the infinite geodesic erosion such as 

= supiElif) e V,g) with E^^if) = /. 
Introduced by Grimaud[9], the geodesic reconstruction is 
called a dynamic filter when the function / is equal to the 
function g plus a constant h: /(*) = g{*) + h. The dynamic 
filter belongs to the category of vertical filter that fills the 
valleys with depth lower than h (see figure 9). 
Introduced by Beucher[4], the geodesic reconstruction is called 
a homotopic transformation when the function / is equal to g 
on the seeds, {sa,...,Sn), and 'oo' on the complementary 
seeds[5] (see figure 10). The homotopic transformation is 
used in the watershed transformation implementation proposed 
by Vincent [18] in order to keep only the most significant 
contours in the areas of interest between the markers. In our 
implementation, the homotopic transformation is done during 
the growing process. 
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Geutlesic dilamiibn of" I" on g 




Image after a dynamie filler 



Fig. 9. The dynamic filter. Before the application of the dynamic filter, there 
are many minima (green bullets). After the application of the dynamic filter, 
there are only two minima. 




the seeds 



Fig. 10. The initial image g with seeds. The function / is equal to ' g{x) 
if X belongs to the seeds and 'oo' if not. 



The classical implementation of the geodisic reconstruction 
is to use directly the formula = sup(i?*(/) Q V,g) 

with Eg{f) = /. Numerically, the recurrence is stopped when 
there is nilpotence, £'*+^(/) = Eg{f). The implementation is 
simple but the complexity is Q{n.k), where n is the number 
of pixels of the image and k is the index of the nilpotence 
condition. 

An alternative to this previous algorithm is an algorithm using 
the SRGPA. The concept of this algorithm is a merging 
procedure. First, a minima procedure is applied on / to extract 
the regional minima (S'i)o<i<g of /. For the convenience, each 
Si is reduced to a single pixel Xi thrown randomly in 5*^. 
The difference with the watershed transformation is that the 
creation of region/ZI is done during the merging procedure. 
At the immersion level equal to level, each region/ZI i is 
created if f{xi) is equal to level and if there is not yet a 
region on the pixel Xi. The last difference is that there is not 
a region boundary to separate two adjacent regions. At every 
growth X of a region, the immersion level is attributed to the 
dynamic function on x, E'^{f){x) = level (see figure 11 and 
algorithm 6). The complexity of this algorithm is 8(71) where 
n is the number of pixels of the image. The application is 
shown on the figure 12. 



I ^ ^ Topographic Eiurf; 



-L-dft- 

Image of dynamique h 




Fig. 11. The process is: image 1) extraction of the regional minima, 
contraction of these regional minima to single pixels Xi, association between 
the chimneys with height f(xi) — g{xi) and the single pixels; image 
2) immersion process: the water enters in the topographic surface by the 
chimneys if there is not a region yet; image 3) catchment basin takes the 
colour of the chimney and at every growth a; of a region, the dynamic image 
takes the immersion level in x; image 4) the red chimney does not create 
region/ZI because the green region is aheady here. Note that there are 3 
minima in the initial image and only two after the dynamic filter 







■ ♦ 

* 

» ■ 













Fig. 12. For the three first serie, the first image is the initial image, the 
second image is the opposite of the distance function of the initial image 
after the application of a dynamic filter. (h=0 for the first serie, h=3 for the 
second, h=10 for the third), the third image is the regional minima of the 
second image, the fourth image is the watershed tranformation restricted by 
the first image on the second image using the minima of the third image like 
a seed. The fourth serie is the same process but in 3D. The first image is 
the initial image. h=0 for the second image, h=3 for the third, h=10 for the 
fourth. 



IV. Conclusion 

In this paper, we implement various algorithms in the 
field of SRGPA. Each implementation is simple and efficient 
using the library Population. When the growing process is 
done at constant velocity with forgetting the past (simulated 
Voronoi tessellation, domain to clusters, regional minima), 
a single queue is sufficient to implement these algorithms. 
When the growing process depends on the topographic surface 
(watershed transformation and dynamic filter) or when an in- 
formation has to be keep during the growing process (distance 
function), the queues number is more than one to implement 
these algorithms. 
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Algorithm 6 Geodesic reconstruction 

Require: /, g V, h //the two images (for the dynamic filter 
f = g + h, the neighbourhood 
II initialization 
Image GR(/ ); 
int level=Q; 

System_Queue s_q( 6{x,i) = ma.x{g{x),level), 
FIFO,.g.max_range() - g.min_range()+l); //n FIFO 
queues. 

Population p (s_q); //create the object Population 

Restricted N ^N; 

Tribe active(V, N); 

iSi)o<t<q= minima(/); 

{xi)o<t<q= rand_pixel((S'i)o<i<g); 

lithe growing process 

. For level = I.min_range() to I.max_range() do 
for i=0 to q do 

// Creation of regionlZI if two conditions 
if ilevel== /(a;i))Apop.X()[a;i].empty()==true then 
r_t= p.growth_tribe(actif ); 
p.growth(a;i, r_t ); 
GR(a;j)=level; 
end if 
end for 

20: s_q.select_queue(ZeweZ); //Select the queue number level 
while s_q.empty()==false do 
(a;,i)=s_q.pop(); 
p.growth(x, i ); 
GR(x)=level; 
end while 
end for 
return GR; 



The application of these algorithms will be present in the two 
next papers of this serie and some new algorithms using the 
SRGPA will be present in a further paper. 

Appendix I 
Definition of distance 

Let be a domain of E, a and b two points of Q.. We call 
geodesic distance b) in A the lowest bound of the length 
of the paths y in Vl, linking a and b. 

Let s be a set. We call the geodesic distance dn{s,b) = 
mm\/aes dx{a,b), the lowest bound of all geodesic distance 
dn{a, b) such as a belongs to s. 
A property of the geodesic distance is: 

dn{a,b) ^ min {dn{a, c) + dnic,b)) 

Vce{a}©v,nn 

The symbol A means the boundary of A. 

We have especially in the discrete space for the norme 1 and 

oo (see figure 13): 

dn{a,b) = min {l + dn{c,b)) 

Vce{a}©Vini2 




Fig. 13. Different paths in the metric .space for the first image and in the 
discrete space for the second image. The red bullets are the neighborhood of 
the point x. The blue bullets are the complementary of the domain f2 



Appendix II 
Summary of the previous articles 

The idea of the first article is to define three objects: Zone 
of Influence (ZI), System of Queues (SQ) and Population. The 
algorithm implementation using SRGPA is focused on the util- 
isation of these three objects. An object ZI is associated to each 
region and localizes a zone on the outer boundary of its region. 
For example, a ZI can be the outer boundary region excluding 
aU other regions. An algorithm using SRGPA is not global 
(no treatment for a block of pixels) but local (the iteration 
is applied pixel by pixel belonging to the ZI). To manage the 
pixel by pixel organisation, a SQ sorts out all pixels belonging 
to ZI depending on the metric and the entering time. It gives 
the possibility to select a pixel following a value of the metric 
and a condition of the entering time. The object population 
links all regions/ZI and permits the (de)growth of regions. 
A pseudo-library, named Population, implements these three 
objects. An algorithm can be implemented easier and faster 
with this Ubrary, fitted for SRGPA. 

The idea of the second article is to give three different growing 
processes, leading up to three different partitions of the space: 

1) one without a boundary region to divide the other 
regions, 

2) another with a boundary region to divide the other 
regions, 

3) the last one does not depend on the seeded region 
initialisation order 
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